Evidence for weak itinerant long-range magnetic correlations in UGe 2 
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Positive muon spin relaxation measurements performed on the ferromagnet UGe2 reveal, in addition 
to the well known localized 5/-electron density responsible for the bulk magnetic properties, the 
existence of itinerant quasi-static magnetic correlations. Their critical dynamics is well described 
by the conventional dipolar Heisenberg model. These correlations involve small magnetic moments. 



PACS numbers : 74.70.Tx, 76.75.+i 

The discovery of superconductivity below 1 K within a 
limited pressure range in the ferromagnet UGe2 [1-4] pro- 
vides an unanticipated example of coexistence of super- 
conductivity and strong ferromagnetism. The electronic 
pairing mechanism needed for superconductivity is be- 
lieved to be magnetic in origin. However, it is amazing 
that ferromagnetically ordered uranium magnetic mo- 
ments with so large magnitude (~ 1.4 /^b at ambient pres- 
sure as deduced from magnetization measurements) are 
directly involved [5]. Since the pairing must involve the 
conduction electrons, it is important to characterize their 
magnetic properties. Because of the restrictions imposed 
by the magnetic form factor, this can not be done by 
diffraction techniques. As the muons localize in intersti- 
tial sites, they have the potentiality to yield information 
on the conduction electrons. Here we show, using the 
muon spin relaxation technique, that UGe2 is actually a 
dual system where two sub-states of / electrons coexist. 
We indeed report the existence at ambient pressure of 
itinerant long-range magnetic correlations with magnetic 
moments of ~ 0.02 /ib and a spectral weight in the mega- 
hertz range. A quantitative understanding of this state 
is moreover reached assuming that these correlations in- 
volve only long wavelength fluctuation modes. 

UGe2 is a ferromagnet with a Curie temperature Tq — 
52 K which crystallizes in the orthorhombic ZrGa2 crystal 
structure (space group Cmmm) [6,7]. Magnetic measure- 
ments indicate a strong magnetocrystalline anisotropy 
[8,9,3] with easy magnetization axis along the a axis. 

We present results obtained by the muon spin relax- 
ation (//SR) technique. Fully polarized muons are im- 
planted into the studied sample. Their spin (1/2) evolves 
in the local magnetic field, Bi OC! until they decay into 
positrons. Since the positron is emitted preferentially in 
the direction of the muon spin at the decay time, it is 
possible to follow the evolution of the muon spin polar- 
ization [10,11]. The measured physical parameter is the 
so-called asymmetry which characterizes the anisotropy 
of the positron emission. Below Tq, if Bi oc has a compo- 



nent perpendicular to the initial muon beam polarization, 
S M (taken parallel to Z), we expect the asymmetry to 
display spontaneous oscillations with an amplitude max- 
imum for Bi oc -L S M . On the other hand, if Bi oc || S^, 
the asymmetry can be written as the product of an ini- 
tial asymmetry related to sample, a s , and the muon spin 
relaxation function, Pz(t), which monitors the dynamics 
of Bi oc . 

UGe2 crystals were grown from a polycrystalline ingot 
using a Czochralski tri-arc technique [8]. We present re- 
sults for two samples. Each consists of pieces cut from 
the crystals, put together to form a disk and glued on a 
silver backing plate. They differ by the orientation (ei- 
ther parallel or perpendicular) of the a axis relative to 
the normal to the sample plane. The measurements were 
performed at the EMU spectrometer of the ISIS facility, 
from 5 K up to 200 K, mostly in zero-field. Additional 
/iSR spectra were recorded with a longitudinal field. 

We found that the temperature dependence of a s for 
S M || a is consistent with Bi oc || a. In agreement 
with that conclusion, a spontaneous muon spin preces- 
sion resulting in wiggles in the asymmetry is observed 
for ± a. Defining Tq as the temperature at which 
the wiggles disappear, we found Tq — 52.49 (2) K. This 
value coincides with the maximum of the relaxation rate 
(to be evidenced below) for S M || a and S M _L a. 

In this letter we focus on the description of data taken 
around the Curie point. 

All the spectra were analyzed as a sum of two compo- 
nents: aP^ p (t) — a s Pz(t) + a Dg . The first component 
describes the /xSR signal from the sample and the second 
accounts for the muons stopped in the background, i.e. 
the cryostat walls and sample holder. In zero-field, for 
all relevant temperatures and for the two orientations of 
Sfj, relative to a, Pz(t) is well described by an exponen- 
tial function: Pz(t) = exp(— Azi) where Xz measures the 
spin-lattice relaxation rate at the muon site. An example 
is shown in Fig. 1. a Dg , which is basically temperature 
independent, was measured for S M _L a and T < Tq as 
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the constant background signal. We got db g = 0.077 [12]. 
For S M || a, it could only be estimated from the sample 
size since the relaxation was never strong enough to mea- 
sure it directly. We took ab g = 0.064. The uncertainty 
on this dbg leads to an uncertainty on the absolute value 
of A Z (S M || a) of - 10%. 
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FIG. 1. Upper panel: examples of piSR spectra recorded 
in zero and longitudinal fields at T = 52.59 (2) K (above 
T c = 52.49 (2) K) for S M _L a. The solid lines are fits assum- 
ing a squared-Lorentzian distribution for the modulus of the 
field at the muon site, Bi oc . The dashed line, which is the 
result of the fit of the zero-field spectrum with an exponential 
relaxation function, can not be distinguished from the solid 
line except above ~ 11 /is. In the lower panel, the comparison 
of the 1.0 mT spectrum with the prediction of an exponential 
fit shows that this model is not valid in longitudinal fields. 
The field dependence at small field of Pz{t) proves that the 
field distribution at the muon site is quasi-static. 

In Fig. 2 we display Xz(T) measured in zero- field for 
S M _L a and S M || a. For both geometries, Xz(T) ex- 
hibits a maximum at Tq. It is due to the critical slowing 
down of the spin dynamics. Surprisingly the anisotropy 
between the orientations is very weak although UGe2 is 
known to be extremely anisotropic [14]. Furthermore we 
show in the following lines that Xz (T) near Tq is quanti- 
tatively understood in the framework of the Heisenberg 
model with dipolar interactions, whereas UGe2 is con- 
sidered as an Ising system. The magnetic signal that 
we observe has therefore a different origin than the well 
documented uranium magnetic state observed e.g. by 
macroscopic measurements. 

Xz{T) has been computed several years ago [15] for the 
critical regime of dipolar Heisenberg ferromagnets and 



has been successfully compared to experiments [15-17]. 
It is based on the derivation of the static and dynamical 
scaling laws from mode coupling theory [18]. The two 
scaling variables at play depend on two material param- 
eters: £o> the magnetic correlation length at T = 2Tc, 
and qn, the dipolar wave vector which is a measure of 
the strength of the exchange interaction relative to the 
dipolar energy. This model initially derived for the para- 
magnetic phase applies also below Tq [16]. 

Specifically, the model predicts that Xz(T) = 
W[a L I L {T) + a T I T (T)} where I L ' T [19] are scaling func- 
tions obtained from mode coupling theory and ol,t are 
parameters determining respectively the amount of longi- 
tudinal (L) and transverse (T) fluctuations probed by the 
measurements. The L,T indices denote the orientation 
relative to the wave vector of the fluctuation mode. ftL.T 
only depend on muon site properties. The result of the fit 
of Xz(T) is shown in the inserts of Fig. 2. The divergence 
of Xz at Tc is strongly reduced by the effect of the dipolar 
interaction [18]. The temperature scale gives the product 
9d£o [15]. For S M _L a, we get g D Co = 0.021 (2), and for 
S„ || a, g D £+ = 0.043 (2) and q D ^ = 0.020 (2). The 
index + (— ) on £o specifies that we consider the param- 
agnetic (ferromagnetic) state. £o(S M || a) > £o(S M T a) 
in the paramagnetic state, suggesting that the magnetic 
correlations are somewhat anisotropic. The fact that 
£o" > £q is an expected feature [10]. The relaxation 
rate scale yields W + a L = 0.140 (4) MHz and W~a L = 
0.20 (2) MHz for S M || a. The transverse contribution 
to X z for both T < T c and T > T c is more difficult to 
estimate since ax is found much lower than a^. Reason- 
able fits are obtained with ot/ol = 0.036 (14). We have 
computed and ax for different possible muon sites and 
found only one site satisfying (Zt < &l/2. This is site 26 
(in Wyckoff notation) of coordinates (0, 1/2, 0) for which 
a L = 1.2486, a T = 0.0386. We then deduce W+ = 0.112 
(3) MHz and W~ = 0.161 (16) MHz. The scale deduced 
from the measurements with S M _L a is about twice as 
large, pointing out again to the weak anisotropy of the 
magnetic correlations. 

In order to further characterize the relaxation near Tc, 
we performed at a given temperature longitudinal field 
measurements for the two orientations of relative to 
a. The field responses for the two geometries are similar. 
An illustration is given in Fig. 1. Surprisingly, the spectra 
are field dependent at extremely low B e xt, proving that 
the probed magnetic fluctuations are quasi-static (fluc- 
tuation rate in the MHz range) and since Xz is small, 
the associated magnetic moment must be small as well. 
Quantitatively, the field dependence of Pz(t) can not be 
described consistently either by a simple exponential re- 
laxation form (see the lower panel of Fig. 1) nor by a 
relaxation function computed with the strong collision 
model assuming an isotropic Gaussian component field 
distribution [20]. On the other hand, the relaxation is 
well explained if we assume that the distribution of the 
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local field at the muon, B\ oc , is squared-Lorentzian [21]. 
We write Pz{t) = Pz{^ot, Vf, t) where AL or character- 
izes the width of the field distribution and Vf its fluctu- 
ation rate [22]. A global fit of the spectra (B cxt = 0, 0.2, 
0.4, 0.6, 0.8, 1.0 and 2.0 mT) taken at a given temper- 
ature is possible. For S M _L a at T = 52.59 (2) K, the 
description of the seven spectra is done with AL or = 70 
fiT and v f = 0.10 MHz. For S M || a at T = 52.47 (2) K 
the two parameters are A Lor = 40 yuT and Vf = 0.50 MHz: 
the zero-field spectra have therefore been recorded in the 
motional narrowing limit (V//(7 M ALor) > 1 where 7^ is 
the muon gyromagnetic ratio; 7 M = 851.6 Mrad s _1 T _1 ). 
This justifies the formalism used to treat Xz(T) close to 
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FIG. 2. 



zero- field for S M _L a and S 



Temperature dependence of Xz measured in 
a in the upper and lower panels, 
respectively. The inserts display Az(T) near Tc- The solid 
and dashed lines are the results of fits for a dipolar Heisenberg 
ferromagnet as explained in the main text. Since Bi oc || a, we 
can not observe a spin-lattice relaxation process for S M _L a in 
the ordered state. The point for S M _L a at (T-T C )/T G = 0.05 
does not fit the critical description, pointing out that it was 
recorded outside the critical region for that geometrical con- 
figuration. The behavior of Az(T) near Tc is typical for a 
ferromagnet in its critical regime. 

We now present an interpretation of our results. 

We first note that the detected fluctuations can not 



arise directly from the localized uranium 5/ electrons 
since Vf would then be in the THz window as estimated 
from Vf ~ k^Tc/h, rather than in the MHz range as mea- 
sured. We also already mentioned that the observed /iSR 
signal has not the properties expected from the known 
macroscopic properties. These apparently conflicting re- 
sults can be understood if the 5/ electrons are viewed 
as two electron subsets. This picture has already been 
argued for UCu 5 [23] and UPd 2 Al 3 [24-28]. However, 
for UGe2 the signatures of both subsets are found at a 
single temperature, the Curie temperature, whereas for 
UCU5 and UPd2Al 3 the temperatures at which the two 
subsets are detected are far apart. So UGe2 presents a 
novel variant of the two electron subset model. Within 
this picture, the anisotropy of the magnetization arises 
from the localized 5/ spectral density and the magnetic 
fluctuations probed by /iSR is a signature of the band-like 
electrons. We do not detect the signature of the localized 
5/ electrons, because of the strong motional narrowing 
of the related relaxation rate. It is observed for the fer- 
romagnet YbNiSn; see Ref. [13]. 

The effect of the dipolar interaction on the quasi-elastic 
linewidth, T(q), of the fluctuations has already been ob- 
served for the weak itinerant ferromagnet M3AI [29]. In 
particular, at criticality T(q) oc q 5 ^ 2 , as expected from 
scaling [18]. Thus it is not completely surprising to de- 
tect its influence on \z{T) for the band-like electrons of 
UGe2. Quantitatively, the data have been described in 
the established framework of critical dynamics [18]. We 
shall now prove that the detected magnetization den- 
sity arises entirely from long wavelength, i.e. small q, 
fluctuations. The magnetic properties of weak itiner- 
ant ferromagnets are explained with the latter hypothesis 
[30,31]. In our model the values of W and Vf and of the 
magnitude of the band- like uranium magnetic moment, 
tou , are controlled by two wavevectors: qu already intro- 
duced and the cut-off wavevector, q c , which sets the up- 
per bound for the wavevector of the fluctuations involved 
in the build-up of the magnetization density. For sim- 
plicity we consider that the magnetic properties of this 
electronic subset are isotropic. We shall detail the anal- 
ysis of the data taken with S M || a. The same approach 
works equally well for the data recorded with _L a. 
As explained below, we get an overall consistent picture 
setting q B = 1.0 x 10~ 3 A" 1 and q c = 0.1 A" 1 . 

The magnetization arising from the conduction elec- 
trons can be viewed as a stochastic variable with a 
variance ((SAi) 2 ). From the fluctuation-dissipation 
(Nyquist's) theorem [32], ((5M) 2 ) obeys the sum rule 
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if the energy of the magnetic fluctuations is smaller 
than the thermal energy. /io is the permeability of 
free space. Assuming an Ornstein-Zernike form for the 
wavevector dependent susceptibility, x(?)i an d since qu 
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is very small, {(SM) 2 ) — 3 k^T q 2 ^ q c /(2ir 2 nn). Since 
m u = v o\/ ((<5-^) 2 ) where Vq is the volume per uranium 
atom (vq = 61.6 A 3 ), we infer m\j — 0.02 at Tq. Inter- 
estingly, the analysis of polarized neutron scattering data 
suggests for the conduction electrons a magnetic moment 
of 0.04 (3) ,ub at low temperature [33]. 

The scale W for Xz can then be computed within the 
framework presented above. Numerically, from Eq. 5.10c 
of Ref. [15], we get W = 0.16 MHz, close to the measured 
values. With the same theory, KT(q) = flq 5 / 2 with ft = 
18 meVA 2 5 at criticality and for small go (see Eq. 4.14b 
of Ref. [15]). Since the measured dynamics is mainly 
driven by the fluctuations at qd [19], we estimate Vj ~ 
r(<?D) = 0.87 MHz, not far from the measured value. 

We now discuss the magnitude of Ai^r • If the distri- 
bution of B\ oc was Gaussian, the zero-field width of the 
distribution would be Acauss = 1.7 mT for muon at site 
2b, and m\j = 0.02 fiQ computed using the Van Vleck 
type formalism of Ref. [20] . However the distribution is 
squared-Lorentzian rather than Gaussian. Such a distri- 
bution is observed in systems with diluted and disordered 
magnetic moments [21]. According to Uemura et al. [22], 
A Lor = ^JVplc/S. 

Gauss where c is the concentration of 
moments at the origin of the distribution. This relation 
leads to c = 1.9 %, consistent with the usual fact that a 
tiny fraction of the total number of valence electrons are 
able to contribute to the magnetic susceptibility. 

From the q^, value we derive the exchange interaction. 
We obtain 2J = fc B T c /4.2 (see Eq. 4.4b of Ref. [15] 
and Ref. [18]). For comparison, the same method gives 
2J/k B T c = 1/11 and 1/20 for metallic Fe and Ni, respec- 
tively. Therefore the evaluation of the exchange energy is 
quite reasonable. From the measured product <Zd£o~, we 
get $ ~ 43 A. This means that the interaction between 
the itinerant magnetic moments is relatively long-range, 
even far outside the critical regime. Although about an 
order of magnitude larger than for conventional ferro- 
magnets, £q~ compares favorably with the neutron result 
for Ni 3 Al: = 24 (9) A [29]. For the same compound, 
we derive from [30] that q c — 0.2 A -1 , a value twice as 
large as found for UGe2 . The moment carried by the itin- 
erant electrons is about four times smaller for UGe2 than 
for N13AI [34]. Nearest neighbor U atoms form zigzag 
chains parallel to a [3]. This may lead to magnetic frus- 
tration and thus explains the disordered nature of the 
distribution of B\ oc . A proper understanding of the ori- 
gin of the squared-Lorentzian distribution requires more 
work. 

One may question the uniqueness of our interpreta- 
tion. We first note that the observed Xz can not arise 
from an impurity phase since the measured critical dy- 
namics occurs right at the well known Curie temperature 
of UGe2 . It could be argued that the observed signal is 
the signature of a weak disorder in the uranium magnetic 
moments. This has already been seen in UAs [35] where 



the ^.SR signal below the Neel temperature has been at- 
tributed to a diluted source of small magnetic moments. 
Their quasi-static nature is related to the absence of spin 
excitations. However, the moments we observe in UGe2 
are quasi-static even above Tq. 

In conclusion we have shown that at ambient pressure 
UGe2 is a dual system where an electronic subset of itin- 
erant states coexist with the subset of localized 5/ elec- 
trons responsible for the well known bulk magnetic prop- 
erties. Its associated magnetic moment is quite small and 
characterized by a very slow spin dynamics. A quanti- 
tative picture for that subset is achieved by assuming 
that only fluctuations at long wavelength are at play. It 
would be of interest to follow the small moment itinerant 
state as a function of pressure to determine whether the 
Cooper's pairs arise from it. However, it seems difficult 
to perform that task with /zSR, unless the spin-lattice 
relaxation rate increases appreciably at high pressure. 
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